3.4.44 \(\int \cos ^5(c+d x) (a+a \sin (c+d x))^m \, dx\) [344]

Optimal. Leaf size=81 \[ \frac {4 (a+a \sin (c+d x))^{3+m}}{a^3 d (3+m)}-\frac {4 (a+a \sin (c+d x))^{4+m}}{a^4 d (4+m)}+\frac {(a+a \sin (c+d x))^{5+m}}{a^5 d (5+m)} \]

[Out]

4*(a+a*sin(d*x+c))^(3+m)/a^3/d/(3+m)-4*(a+a*sin(d*x+c))^(4+m)/a^4/d/(4+m)+(a+a*sin(d*x+c))^(5+m)/a^5/d/(5+m)

________________________________________________________________________________________

Rubi [A]
time = 0.05, antiderivative size = 81, normalized size of antiderivative = 1.00, number of steps used = 3, number of rules used = 2, integrand size = 21, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.095, Rules used = {2746, 45} \begin {gather*} \frac {(a \sin (c+d x)+a)^{m+5}}{a^5 d (m+5)}-\frac {4 (a \sin (c+d x)+a)^{m+4}}{a^4 d (m+4)}+\frac {4 (a \sin (c+d x)+a)^{m+3}}{a^3 d (m+3)} \end {gather*}

Antiderivative was successfully verified.

[In]

Int[Cos[c + d*x]^5*(a + a*Sin[c + d*x])^m,x]

[Out]

(4*(a + a*Sin[c + d*x])^(3 + m))/(a^3*d*(3 + m)) - (4*(a + a*Sin[c + d*x])^(4 + m))/(a^4*d*(4 + m)) + (a + a*S
in[c + d*x])^(5 + m)/(a^5*d*(5 + m))

Rule 45

Int[((a_.) + (b_.)*(x_))^(m_.)*((c_.) + (d_.)*(x_))^(n_.), x_Symbol] :> Int[ExpandIntegrand[(a + b*x)^m*(c + d
*x)^n, x], x] /; FreeQ[{a, b, c, d, n}, x] && NeQ[b*c - a*d, 0] && IGtQ[m, 0] && ( !IntegerQ[n] || (EqQ[c, 0]
&& LeQ[7*m + 4*n + 4, 0]) || LtQ[9*m + 5*(n + 1), 0] || GtQ[m + n + 2, 0])

Rule 2746

Int[cos[(e_.) + (f_.)*(x_)]^(p_.)*((a_) + (b_.)*sin[(e_.) + (f_.)*(x_)])^(m_.), x_Symbol] :> Dist[1/(b^p*f), S
ubst[Int[(a + x)^(m + (p - 1)/2)*(a - x)^((p - 1)/2), x], x, b*Sin[e + f*x]], x] /; FreeQ[{a, b, e, f, m}, x]
&& IntegerQ[(p - 1)/2] && EqQ[a^2 - b^2, 0] && (GeQ[p, -1] ||  !IntegerQ[m + 1/2])

Rubi steps

\begin {align*} \int \cos ^5(c+d x) (a+a \sin (c+d x))^m \, dx &=\frac {\text {Subst}\left (\int (a-x)^2 (a+x)^{2+m} \, dx,x,a \sin (c+d x)\right )}{a^5 d}\\ &=\frac {\text {Subst}\left (\int \left (4 a^2 (a+x)^{2+m}-4 a (a+x)^{3+m}+(a+x)^{4+m}\right ) \, dx,x,a \sin (c+d x)\right )}{a^5 d}\\ &=\frac {4 (a+a \sin (c+d x))^{3+m}}{a^3 d (3+m)}-\frac {4 (a+a \sin (c+d x))^{4+m}}{a^4 d (4+m)}+\frac {(a+a \sin (c+d x))^{5+m}}{a^5 d (5+m)}\\ \end {align*}

________________________________________________________________________________________

Mathematica [A]
time = 1.48, size = 100, normalized size = 1.23 \begin {gather*} -\frac {\left (\cos \left (\frac {1}{2} (c+d x)\right )+\sin \left (\frac {1}{2} (c+d x)\right )\right )^6 (a (1+\sin (c+d x)))^m \left (-76-29 m-3 m^2+\left (12+7 m+m^2\right ) \cos (2 (c+d x))+4 \left (18+9 m+m^2\right ) \sin (c+d x)\right )}{2 d (3+m) (4+m) (5+m)} \end {gather*}

Antiderivative was successfully verified.

[In]

Integrate[Cos[c + d*x]^5*(a + a*Sin[c + d*x])^m,x]

[Out]

-1/2*((Cos[(c + d*x)/2] + Sin[(c + d*x)/2])^6*(a*(1 + Sin[c + d*x]))^m*(-76 - 29*m - 3*m^2 + (12 + 7*m + m^2)*
Cos[2*(c + d*x)] + 4*(18 + 9*m + m^2)*Sin[c + d*x]))/(d*(3 + m)*(4 + m)*(5 + m))

________________________________________________________________________________________

Maple [F]
time = 0.30, size = 0, normalized size = 0.00 \[\int \left (\cos ^{5}\left (d x +c \right )\right ) \left (a +a \sin \left (d x +c \right )\right )^{m}\, dx\]

Verification of antiderivative is not currently implemented for this CAS.

[In]

int(cos(d*x+c)^5*(a+a*sin(d*x+c))^m,x)

[Out]

int(cos(d*x+c)^5*(a+a*sin(d*x+c))^m,x)

________________________________________________________________________________________

Maxima [B] Leaf count of result is larger than twice the leaf count of optimal. 266 vs. \(2 (81) = 162\).
time = 0.28, size = 266, normalized size = 3.28 \begin {gather*} \frac {\frac {{\left ({\left (m^{4} + 10 \, m^{3} + 35 \, m^{2} + 50 \, m + 24\right )} a^{m} \sin \left (d x + c\right )^{5} + {\left (m^{4} + 6 \, m^{3} + 11 \, m^{2} + 6 \, m\right )} a^{m} \sin \left (d x + c\right )^{4} - 4 \, {\left (m^{3} + 3 \, m^{2} + 2 \, m\right )} a^{m} \sin \left (d x + c\right )^{3} + 12 \, {\left (m^{2} + m\right )} a^{m} \sin \left (d x + c\right )^{2} - 24 \, a^{m} m \sin \left (d x + c\right ) + 24 \, a^{m}\right )} {\left (\sin \left (d x + c\right ) + 1\right )}^{m}}{m^{5} + 15 \, m^{4} + 85 \, m^{3} + 225 \, m^{2} + 274 \, m + 120} - \frac {2 \, {\left ({\left (m^{2} + 3 \, m + 2\right )} a^{m} \sin \left (d x + c\right )^{3} + {\left (m^{2} + m\right )} a^{m} \sin \left (d x + c\right )^{2} - 2 \, a^{m} m \sin \left (d x + c\right ) + 2 \, a^{m}\right )} {\left (\sin \left (d x + c\right ) + 1\right )}^{m}}{m^{3} + 6 \, m^{2} + 11 \, m + 6} + \frac {{\left (a \sin \left (d x + c\right ) + a\right )}^{m + 1}}{a {\left (m + 1\right )}}}{d} \end {gather*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(cos(d*x+c)^5*(a+a*sin(d*x+c))^m,x, algorithm="maxima")

[Out]

(((m^4 + 10*m^3 + 35*m^2 + 50*m + 24)*a^m*sin(d*x + c)^5 + (m^4 + 6*m^3 + 11*m^2 + 6*m)*a^m*sin(d*x + c)^4 - 4
*(m^3 + 3*m^2 + 2*m)*a^m*sin(d*x + c)^3 + 12*(m^2 + m)*a^m*sin(d*x + c)^2 - 24*a^m*m*sin(d*x + c) + 24*a^m)*(s
in(d*x + c) + 1)^m/(m^5 + 15*m^4 + 85*m^3 + 225*m^2 + 274*m + 120) - 2*((m^2 + 3*m + 2)*a^m*sin(d*x + c)^3 + (
m^2 + m)*a^m*sin(d*x + c)^2 - 2*a^m*m*sin(d*x + c) + 2*a^m)*(sin(d*x + c) + 1)^m/(m^3 + 6*m^2 + 11*m + 6) + (a
*sin(d*x + c) + a)^(m + 1)/(a*(m + 1)))/d

________________________________________________________________________________________

Fricas [A]
time = 0.38, size = 102, normalized size = 1.26 \begin {gather*} \frac {{\left ({\left (m^{2} + 3 \, m\right )} \cos \left (d x + c\right )^{4} + 8 \, m \cos \left (d x + c\right )^{2} + {\left ({\left (m^{2} + 7 \, m + 12\right )} \cos \left (d x + c\right )^{4} + 8 \, {\left (m + 2\right )} \cos \left (d x + c\right )^{2} + 32\right )} \sin \left (d x + c\right ) + 32\right )} {\left (a \sin \left (d x + c\right ) + a\right )}^{m}}{d m^{3} + 12 \, d m^{2} + 47 \, d m + 60 \, d} \end {gather*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(cos(d*x+c)^5*(a+a*sin(d*x+c))^m,x, algorithm="fricas")

[Out]

((m^2 + 3*m)*cos(d*x + c)^4 + 8*m*cos(d*x + c)^2 + ((m^2 + 7*m + 12)*cos(d*x + c)^4 + 8*(m + 2)*cos(d*x + c)^2
+ 32)*sin(d*x + c) + 32)*(a*sin(d*x + c) + a)^m/(d*m^3 + 12*d*m^2 + 47*d*m + 60*d)

________________________________________________________________________________________

Sympy [B] Leaf count of result is larger than twice the leaf count of optimal. 5534 vs. \(2 (68) = 136\).
time = 42.01, size = 5534, normalized size = 68.32 \begin {gather*} \text {Too large to display} \end {gather*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(cos(d*x+c)**5*(a+a*sin(d*x+c))**m,x)

[Out]

Piecewise((x*(a*sin(c) + a)**m*cos(c)**5, Eq(d, 0)), (12*log(sin(c + d*x) + 1)*sin(c + d*x)**4/(12*a**5*d*sin(
c + d*x)**4 + 48*a**5*d*sin(c + d*x)**3 + 72*a**5*d*sin(c + d*x)**2 + 48*a**5*d*sin(c + d*x) + 12*a**5*d) + 48
*log(sin(c + d*x) + 1)*sin(c + d*x)**3/(12*a**5*d*sin(c + d*x)**4 + 48*a**5*d*sin(c + d*x)**3 + 72*a**5*d*sin(
c + d*x)**2 + 48*a**5*d*sin(c + d*x) + 12*a**5*d) + 72*log(sin(c + d*x) + 1)*sin(c + d*x)**2/(12*a**5*d*sin(c
+ d*x)**4 + 48*a**5*d*sin(c + d*x)**3 + 72*a**5*d*sin(c + d*x)**2 + 48*a**5*d*sin(c + d*x) + 12*a**5*d) + 48*l
og(sin(c + d*x) + 1)*sin(c + d*x)/(12*a**5*d*sin(c + d*x)**4 + 48*a**5*d*sin(c + d*x)**3 + 72*a**5*d*sin(c + d
*x)**2 + 48*a**5*d*sin(c + d*x) + 12*a**5*d) + 12*log(sin(c + d*x) + 1)/(12*a**5*d*sin(c + d*x)**4 + 48*a**5*d
*sin(c + d*x)**3 + 72*a**5*d*sin(c + d*x)**2 + 48*a**5*d*sin(c + d*x) + 12*a**5*d) + 20*sin(c + d*x)**3/(12*a*
*5*d*sin(c + d*x)**4 + 48*a**5*d*sin(c + d*x)**3 + 72*a**5*d*sin(c + d*x)**2 + 48*a**5*d*sin(c + d*x) + 12*a**
5*d) + 6*sin(c + d*x)**2*cos(c + d*x)**2/(12*a**5*d*sin(c + d*x)**4 + 48*a**5*d*sin(c + d*x)**3 + 72*a**5*d*si
n(c + d*x)**2 + 48*a**5*d*sin(c + d*x) + 12*a**5*d) + 56*sin(c + d*x)**2/(12*a**5*d*sin(c + d*x)**4 + 48*a**5*
d*sin(c + d*x)**3 + 72*a**5*d*sin(c + d*x)**2 + 48*a**5*d*sin(c + d*x) + 12*a**5*d) + 8*sin(c + d*x)*cos(c + d
*x)**2/(12*a**5*d*sin(c + d*x)**4 + 48*a**5*d*sin(c + d*x)**3 + 72*a**5*d*sin(c + d*x)**2 + 48*a**5*d*sin(c +
d*x) + 12*a**5*d) + 52*sin(c + d*x)/(12*a**5*d*sin(c + d*x)**4 + 48*a**5*d*sin(c + d*x)**3 + 72*a**5*d*sin(c +
d*x)**2 + 48*a**5*d*sin(c + d*x) + 12*a**5*d) - 3*cos(c + d*x)**4/(12*a**5*d*sin(c + d*x)**4 + 48*a**5*d*sin(
c + d*x)**3 + 72*a**5*d*sin(c + d*x)**2 + 48*a**5*d*sin(c + d*x) + 12*a**5*d) + 2*cos(c + d*x)**2/(12*a**5*d*s
in(c + d*x)**4 + 48*a**5*d*sin(c + d*x)**3 + 72*a**5*d*sin(c + d*x)**2 + 48*a**5*d*sin(c + d*x) + 12*a**5*d) +
16/(12*a**5*d*sin(c + d*x)**4 + 48*a**5*d*sin(c + d*x)**3 + 72*a**5*d*sin(c + d*x)**2 + 48*a**5*d*sin(c + d*x
) + 12*a**5*d), Eq(m, -5)), (-12*log(sin(c + d*x) + 1)*sin(c + d*x)**3/(3*a**4*d*sin(c + d*x)**3 + 9*a**4*d*si
n(c + d*x)**2 + 9*a**4*d*sin(c + d*x) + 3*a**4*d) - 36*log(sin(c + d*x) + 1)*sin(c + d*x)**2/(3*a**4*d*sin(c +
d*x)**3 + 9*a**4*d*sin(c + d*x)**2 + 9*a**4*d*sin(c + d*x) + 3*a**4*d) - 36*log(sin(c + d*x) + 1)*sin(c + d*x
)/(3*a**4*d*sin(c + d*x)**3 + 9*a**4*d*sin(c + d*x)**2 + 9*a**4*d*sin(c + d*x) + 3*a**4*d) - 12*log(sin(c + d*
x) + 1)/(3*a**4*d*sin(c + d*x)**3 + 9*a**4*d*sin(c + d*x)**2 + 9*a**4*d*sin(c + d*x) + 3*a**4*d) + 8*sin(c + d
*x)**4/(3*a**4*d*sin(c + d*x)**3 + 9*a**4*d*sin(c + d*x)**2 + 9*a**4*d*sin(c + d*x) + 3*a**4*d) + 4*sin(c + d*
x)**2*cos(c + d*x)**2/(3*a**4*d*sin(c + d*x)**3 + 9*a**4*d*sin(c + d*x)**2 + 9*a**4*d*sin(c + d*x) + 3*a**4*d)
- 52*sin(c + d*x)**2/(3*a**4*d*sin(c + d*x)**3 + 9*a**4*d*sin(c + d*x)**2 + 9*a**4*d*sin(c + d*x) + 3*a**4*d)
+ 6*sin(c + d*x)*cos(c + d*x)**2/(3*a**4*d*sin(c + d*x)**3 + 9*a**4*d*sin(c + d*x)**2 + 9*a**4*d*sin(c + d*x)
+ 3*a**4*d) - 72*sin(c + d*x)/(3*a**4*d*sin(c + d*x)**3 + 9*a**4*d*sin(c + d*x)**2 + 9*a**4*d*sin(c + d*x) +
3*a**4*d) - cos(c + d*x)**4/(3*a**4*d*sin(c + d*x)**3 + 9*a**4*d*sin(c + d*x)**2 + 9*a**4*d*sin(c + d*x) + 3*a
**4*d) + 2*cos(c + d*x)**2/(3*a**4*d*sin(c + d*x)**3 + 9*a**4*d*sin(c + d*x)**2 + 9*a**4*d*sin(c + d*x) + 3*a*
*4*d) - 28/(3*a**4*d*sin(c + d*x)**3 + 9*a**4*d*sin(c + d*x)**2 + 9*a**4*d*sin(c + d*x) + 3*a**4*d), Eq(m, -4)
), (8*log(tan(c/2 + d*x/2) + 1)*tan(c/2 + d*x/2)**4/(a**3*d*tan(c/2 + d*x/2)**4 + 2*a**3*d*tan(c/2 + d*x/2)**2
+ a**3*d) + 16*log(tan(c/2 + d*x/2) + 1)*tan(c/2 + d*x/2)**2/(a**3*d*tan(c/2 + d*x/2)**4 + 2*a**3*d*tan(c/2 +
d*x/2)**2 + a**3*d) + 8*log(tan(c/2 + d*x/2) + 1)/(a**3*d*tan(c/2 + d*x/2)**4 + 2*a**3*d*tan(c/2 + d*x/2)**2
+ a**3*d) - 4*log(tan(c/2 + d*x/2)**2 + 1)*tan(c/2 + d*x/2)**4/(a**3*d*tan(c/2 + d*x/2)**4 + 2*a**3*d*tan(c/2
+ d*x/2)**2 + a**3*d) - 8*log(tan(c/2 + d*x/2)**2 + 1)*tan(c/2 + d*x/2)**2/(a**3*d*tan(c/2 + d*x/2)**4 + 2*a**
3*d*tan(c/2 + d*x/2)**2 + a**3*d) - 4*log(tan(c/2 + d*x/2)**2 + 1)/(a**3*d*tan(c/2 + d*x/2)**4 + 2*a**3*d*tan(
c/2 + d*x/2)**2 + a**3*d) - 6*tan(c/2 + d*x/2)**3/(a**3*d*tan(c/2 + d*x/2)**4 + 2*a**3*d*tan(c/2 + d*x/2)**2 +
a**3*d) + 2*tan(c/2 + d*x/2)**2/(a**3*d*tan(c/2 + d*x/2)**4 + 2*a**3*d*tan(c/2 + d*x/2)**2 + a**3*d) - 6*tan(
c/2 + d*x/2)/(a**3*d*tan(c/2 + d*x/2)**4 + 2*a**3*d*tan(c/2 + d*x/2)**2 + a**3*d), Eq(m, -3)), (6*tan(c/2 + d*
x/2)**5/(3*a**2*d*tan(c/2 + d*x/2)**6 + 9*a**2*d*tan(c/2 + d*x/2)**4 + 9*a**2*d*tan(c/2 + d*x/2)**2 + 3*a**2*d
) - 12*tan(c/2 + d*x/2)**4/(3*a**2*d*tan(c/2 + d*x/2)**6 + 9*a**2*d*tan(c/2 + d*x/2)**4 + 9*a**2*d*tan(c/2 + d
*x/2)**2 + 3*a**2*d) + 20*tan(c/2 + d*x/2)**3/(3*a**2*d*tan(c/2 + d*x/2)**6 + 9*a**2*d*tan(c/2 + d*x/2)**4 + 9
*a**2*d*tan(c/2 + d*x/2)**2 + 3*a**2*d) - 12*tan(c/2 + d*x/2)**2/(3*a**2*d*tan(c/2 + d*x/2)**6 + 9*a**2*d*tan(
c/2 + d*x/2)**4 + 9*a**2*d*tan(c/2 + d*x/2)**2 + 3*a**2*d) + 6*tan(c/2 + d*x/2)/(3*a**2*d*tan(c/2 + d*x/2)**6
+ 9*a**2*d*tan(c/2 + d*x/2)**4 + 9*a**2*d*tan(c...

________________________________________________________________________________________

Giac [B] Leaf count of result is larger than twice the leaf count of optimal. 294 vs. \(2 (81) = 162\).
time = 4.38, size = 294, normalized size = 3.63 \begin {gather*} \frac {{\left (a \sin \left (d x + c\right ) + a\right )}^{5} {\left (a \sin \left (d x + c\right ) + a\right )}^{m} m^{2} - 4 \, {\left (a \sin \left (d x + c\right ) + a\right )}^{4} {\left (a \sin \left (d x + c\right ) + a\right )}^{m} a m^{2} + 4 \, {\left (a \sin \left (d x + c\right ) + a\right )}^{3} {\left (a \sin \left (d x + c\right ) + a\right )}^{m} a^{2} m^{2} + 7 \, {\left (a \sin \left (d x + c\right ) + a\right )}^{5} {\left (a \sin \left (d x + c\right ) + a\right )}^{m} m - 32 \, {\left (a \sin \left (d x + c\right ) + a\right )}^{4} {\left (a \sin \left (d x + c\right ) + a\right )}^{m} a m + 36 \, {\left (a \sin \left (d x + c\right ) + a\right )}^{3} {\left (a \sin \left (d x + c\right ) + a\right )}^{m} a^{2} m + 12 \, {\left (a \sin \left (d x + c\right ) + a\right )}^{5} {\left (a \sin \left (d x + c\right ) + a\right )}^{m} - 60 \, {\left (a \sin \left (d x + c\right ) + a\right )}^{4} {\left (a \sin \left (d x + c\right ) + a\right )}^{m} a + 80 \, {\left (a \sin \left (d x + c\right ) + a\right )}^{3} {\left (a \sin \left (d x + c\right ) + a\right )}^{m} a^{2}}{{\left (a^{4} m^{3} + 12 \, a^{4} m^{2} + 47 \, a^{4} m + 60 \, a^{4}\right )} a d} \end {gather*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(cos(d*x+c)^5*(a+a*sin(d*x+c))^m,x, algorithm="giac")

[Out]

((a*sin(d*x + c) + a)^5*(a*sin(d*x + c) + a)^m*m^2 - 4*(a*sin(d*x + c) + a)^4*(a*sin(d*x + c) + a)^m*a*m^2 + 4
*(a*sin(d*x + c) + a)^3*(a*sin(d*x + c) + a)^m*a^2*m^2 + 7*(a*sin(d*x + c) + a)^5*(a*sin(d*x + c) + a)^m*m - 3
2*(a*sin(d*x + c) + a)^4*(a*sin(d*x + c) + a)^m*a*m + 36*(a*sin(d*x + c) + a)^3*(a*sin(d*x + c) + a)^m*a^2*m +
12*(a*sin(d*x + c) + a)^5*(a*sin(d*x + c) + a)^m - 60*(a*sin(d*x + c) + a)^4*(a*sin(d*x + c) + a)^m*a + 80*(a
*sin(d*x + c) + a)^3*(a*sin(d*x + c) + a)^m*a^2)/((a^4*m^3 + 12*a^4*m^2 + 47*a^4*m + 60*a^4)*a*d)

________________________________________________________________________________________

Mupad [B]
time = 1.97, size = 195, normalized size = 2.41 \begin {gather*} \frac {{\left (a\,\left (\sin \left (c+d\,x\right )+1\right )\right )}^m\,\left (82\,m+600\,\sin \left (c+d\,x\right )+100\,\sin \left (3\,c+3\,d\,x\right )+12\,\sin \left (5\,c+5\,d\,x\right )+46\,m\,\sin \left (c+d\,x\right )+88\,m\,\cos \left (2\,c+2\,d\,x\right )+6\,m\,\cos \left (4\,c+4\,d\,x\right )+53\,m\,\sin \left (3\,c+3\,d\,x\right )+7\,m\,\sin \left (5\,c+5\,d\,x\right )+2\,m^2\,\sin \left (c+d\,x\right )+6\,m^2+8\,m^2\,\cos \left (2\,c+2\,d\,x\right )+2\,m^2\,\cos \left (4\,c+4\,d\,x\right )+3\,m^2\,\sin \left (3\,c+3\,d\,x\right )+m^2\,\sin \left (5\,c+5\,d\,x\right )+512\right )}{16\,d\,\left (m^3+12\,m^2+47\,m+60\right )} \end {gather*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

int(cos(c + d*x)^5*(a + a*sin(c + d*x))^m,x)

[Out]

((a*(sin(c + d*x) + 1))^m*(82*m + 600*sin(c + d*x) + 100*sin(3*c + 3*d*x) + 12*sin(5*c + 5*d*x) + 46*m*sin(c +
d*x) + 88*m*cos(2*c + 2*d*x) + 6*m*cos(4*c + 4*d*x) + 53*m*sin(3*c + 3*d*x) + 7*m*sin(5*c + 5*d*x) + 2*m^2*si
n(c + d*x) + 6*m^2 + 8*m^2*cos(2*c + 2*d*x) + 2*m^2*cos(4*c + 4*d*x) + 3*m^2*sin(3*c + 3*d*x) + m^2*sin(5*c +
5*d*x) + 512))/(16*d*(47*m + 12*m^2 + m^3 + 60))

________________________________________________________________________________________